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ABSTRACT 

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation 
and study its internal properties to the present epoch. We use a set of smoothed particle hydrodynamic (SPH) 
simulations with a high dynamical range and force resolution that include cooling, star formation, supernovae 
(SN) feedback and a redshift dependent UV background. We compare results from a ACDM simulation with a 
AWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force 
resolution in both the gas and dark matter components play an important role in solving the angular momentum 
catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. Hence, 
a large disk forms without the need of strong energy injection, the z = galaxies lie close to the I-band Tully- 
Fisher relation, and the stellar material in the disk component has a final specific angular momentum equal 
to 40% and 90% of that of the dark halo in the ACDM and AWDM models respectively. If rescaled to the 
Milky Way, the ACDM galaxy has an overabundance of satellites, with a total mass in the stellar halo 40% that 
in the bulge+disk system. The AWDM galaxy has a drastically reduced satellite population and a negligible 
stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites 
possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around 
the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds 
up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating 
a realistic stellar disk in a massive galaxy within the ACDM scenario, energy injection emerges as necessary 
ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. 
Subject headings: galaxies: evolution - galaxies: formation - methods: numerical 



1. INTRODUCTION 

Within the current paradigm of galaxy formation, realistic 
disks form only if gas retains most of its angular momentum 
gained by torques from nearby structures while it cools at the 
center of cold dark matter (CDM) halos (White & Rees 1978, 
Fall & Efstathiou 1980, Fall 1983, Mo, Mao & White, 1998). 
First simulations of galaxy formation that included star for- 
mation (Lake & Carlberg 1988, Katz 1992) provided strong 
evidence that hierarchical models do create rotationally sup- 
ported stellar systems. However, simulations of galaxy for- 
mation in a full cosmological context have not yet been able 
to form realistic disk galaxies: dynamical friction suffered 
by dense gaseous lumps and subsequent catastrophic angu- 
lar momentum loss caused typical disk scale lengths to come 
short of those observed (Navarro & White 1994). In addi- 
tion, high resolution simulations of dark matter halos in CDM 
models have far more substructure ("satellites") than observa- 
tions; encounters with these satellites could destroy the stellar 
disk component of spiral galaxies (Moore et al. 1999). 

If these dark subhalos were the main problem, any mech- 
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anism capable of reducing the lumpiness of galaxy assem- 
bly would lead to realistic disks. However resolution sim- 
ulations (N ~ 10 4 ) also suffer from excessive artificial vis- 
cosity that would overestimate the angular momentum trans- 
fer, as was early recognized in previous works (e.g (Sommer- 
Larsen, Gelato & Vedel 1999). Poorly resolved disk will also 
suffer from two-body effects which will lead to substantial 
disk angular momentum loss even after its eventual formation 
(Mayer et al. 2001a). 

Together, the failure to build galactic disks similar to those 
observed and the overabundance of substructure posed a 
formidable challenge to the ACDM paradigm, suggesting that 
a fundamental ingredient is missing from our understand- 
ing of the formation of galactic and subgalactic structures. 
As possible solutions, strong SN feedback (e.g Thacker & 
Couchman 2000 (TC) and 2001 (TC01)), an external UV 
background produced by QSOs and massive stars (Quinn, 
Katz, & Efstathiou 1996, Benson et al. 2002, Somerville 
2002) and alternatives to ACDM (Spergel & Steinhardt 2000, 
among others) have been advocated. Among alternative dark 
matter models, Warm Dark Matter (WDM) (Pagels & Pri- 
mack 1982) emerged as a good candidate as its lack of power 
at small scales could prove beneficial in reducing the num- 
ber of subhalos and allowing a smoother accretion of gas, 
forming larger disks (Bode, Ostriker & Turok 2001, Colin, 
Avila-Reese, & Valenzuela 2000). Some of these solutions 
have been explored in recent studies (e.g. TC01, Sommer- 
Larsen & Dolgov 2001, Sommer-Larsen, Gotz & Portinari 
2003, Abadi et al 2003). 

These most recent works have been successful in forming 
a rotationally supported stellar component although in many 
cases a massive stellar spheroid is formed as well. In runs 
that use less than a few thousand particles to represent the 
virialized part of a galaxy dark matter halo (e.g Eke, Navarro 
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& Steinmetz 2001) the rotating stellar component does not al- 
low to distinguish a disk. In TC01, where ~ 1.5 10 4 dark mat- 
ter (DM) particles with the halo virial radius where used, the 
"bulge" component is still almost twice as massive as the thin 
disk. The bulge component decreases to only 50% of the stel- 
lar disk in the simulation presented by Abadi et al (2003) that 
had 3.6 10 4 DM particles within the virial radius of their sim- 
ulated galaxy. While these results are encouraging, previous 
simulations do not clarify if stronger feedback was necessary, 
or just sufficient to avoid systematic loss of angular momen- 
tum during the formation of galactic disks. Runs with larger 
particle numbers also had larger disks, suggesting that reso- 
lution effects still play a major role. Important differences in 
the results might also arise from the different merging histo- 
ries of each individual halo simulated, the spin of their halos 
(e.g TC01 simulated a halo with a larger than average spin 
which in turn makes it easier to form larger disks) and the 
different SN feedback recipes used in various works. 

We have performed a new set of smoothed particle hydrody- 
namic (SPH) simulations that use a high number of particles 
and high force resolution to follow the formation and internal 
structure of the stellar component of galaxies. These simula- 
tions implement cooling, star formation, supernova feedback 
and a UV background. In this paper, we describe the for- 
mation and evolution of a realistic disk dominated galaxies 
formed in a cosmological numerical simulation and compare 
a high-resolution spiral galaxy formed in 2 different (ACDM 
and AWDM) flat cosmologies. Rather than simulate a large 
sample of galaxies, our strategy is to achieve the highest pos- 
sible resolution with the available computing resources to 
avoid unwanted numerical systematics. Our aim is to find a 
good, clear cut case where a realistic disk forms in a fully cos- 
mological environment rather than exploring the full range of 
galactic morphologies. We also evaluate the various numer- 
ical and physical processes responsible for the angular mo- 
mentum loss of the gaseous and stellar components. We de- 
scribe the initial conditions in section 2, we present our results 
and discuss their implications in sections 3 and 4 respectively. 

2. SIMULATION PARAMETERS AND INITIAL CONDITIONS 

We adopted a flat A-dominated cosmology: fio = 0.3, 
A=0.7, h = 0.7, erg = 1, shape parameter T = 0.21, and flt> = 
0.039 (Perlmutter et al. 1997, Efstathiou et al. 2002). Power 
spectra were calculated using the CMBFAST code to generate 
transfer functions (Zaldarriaga & Seljak 2000). The AWDM 
model can be characterized in terms of the reduction in fluc- 
tuations on dwarf galaxy scales. For our choice of 2 keV neu- 
trinos the fluctuations on scales smaller than <~ 10 9 5 M Q are 
vastly reduced. We retained the phases of the waves that de- 
scribe the initial density field. These allowed us to discrim- 
inate between the different effects of the two power spectra 
used. 

From a 100 Mpc box, low resolution simulation, we se- 
lected and resimulated (Katz & White 1993) at higher res- 
olution a "typical" spiral galaxy candidate halo with mass 3 
x 10 12 M Q and a quiet merger history after z ~ 2.5. The halo is 
relatively isolated and its formation time (defined as the main 
progenitor achieving 50% of its final mass) is z = 0.75. Its 
spin parameter, A is 0.035 (0.03 for A WDM) close to the 
average value for cosmic halos (Gardner 2001). Note that 
due to the finite size of the box the actual erg in the simula- 
tion volume is actually ~ 0.9 at the present time. The above 
constraints are consistent with the Milky Way forming sev- 
eral Gyrs ago (Wyse 2002) and predictions from semianalyt- 
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FIG. 1. — Stellar surface density profiles for the different stellar compo- 
nents at z = for the ACDM galaxy. Filled dots are for the total stellar surface 
density profile, open dots for the spheroid only (stars older than 10 Gyr). The 
black lines are exponential fit to the disk component, with scale length 3 kpc 
(3.6 kpc for AWDM). The short dashed lines are exponential fits to the cen- 
tral stellar mass distribution, with scale lengths of 1 kpc for both galaxies 
(the bar and the bulge are both contributing), and the long dashed line is a de 
Vaucouleurs fit to the spheroid, with scale lengths of 800 pc ( ACDM galaxy) 
and 500 pc (AWDM galaxy) 
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FIG. 2. — Stellar surface density profiles for the different stellar components 
at z = for the AWDM galaxy. Filled dots are for the total stellar surface 
density profile, open dots for the spheroid only (stars older than 10 Gyr). The 
black lines are exponential fit to the disk component, with scale length 3.6 
kpc. The short dashed lines are exponential fits to the central stellar mass 
distribution, with scale lengths of 1 kpc for both galaxies (the bar and the 
bulge are both contributing), and the long dashed line is a de Vaucouleurs fit 
to the spheroid, with scale lengths of 500 pc 
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ical models (Baugh, Cole & Frenk 1996). A large starting 
box is important as lack of large scale power would signifi- 
cantly reduce the amount of torques on collapsing halos. To 
run the simulations, we used a parallel treecode+SPH with 
multiple timesteps: GASOLINE described in detail in Wads- 
ley, Stadel & Quinn (2003) and Stadel, Wadsley & Richard- 
son (2001). The code treats artificial viscosity as suggested 
in Balsara (1997). The energy equation is computed asym- 
metrically (Springel & Hernquist 2002, Evrard 1988, Mon- 
aghan 1992 ). This approach avoids the energy conservation 
and negative energy problems of the arithmetic and geometric 
forms and converges to the high resolution answer faster than 
other proposed methods (Benz 1990, Springel and Hernquist 
2002). 

The high resolution region has dark matter and gas parti- 
cles of mass 2.32 x 10 7 and 3.44 x 10 6 M© respectively. The 
rest of the simulation box was re-sampled at increasing parti- 
cle masses for a total of nearly 1 . 1 million dark matter parti- 
cles and 2.2x 10 5 gas particles. For all particles species in the 
high resolution region, the gravitational spline softening, e(z), 
evolved comovingly starting at z=100 until z=9, and remained 
fixed at 1 kpc from z=9 to the present. This choice of e(z) re- 
duces two body relaxation at high z in small halos. The value 
chosen is a good compromise between reducing two body re- 
laxation and ensuring that disk scale lengths and the central 
part of dark matter halos will be spatially resolved (see Die- 
mand et al. 2002 for a number of relevant tests). Moreover, a 
constant softening in physical coordinates is more appropriate 
at low redshift when stellar structures detached from the Hub- 
ble flow. Integration parameter values were chosen as sug- 
gested in Moore et al. (1998) and then confirmed in Power et 
al. (2002). Most of the stars formed in the simulation required 
32,768 integration steps to z = 0, with a significant fraction of 
them requiring four times as many steps. At comparable or 
better force resolution, particle number within the virial radius 
is three times as in Navarro & Steinmetz 2002, eight times as 
in TC01, three times as in Abadi et al. 2002 and compara- 
ble to the highest resolution run reported in Sommer-Larsen, 
Gotz & Portinari 2003). 

Two main SPH simulations were run, one for each cosmol- 
ogy. Both runs included: Compton and radiative cooling, as- 
suming a gas of primordial composition; star formation and 
SN (type I&II) feedback, treated following the prescription 
described by Katz (1992), where stars spawn from cold, Jeans 
unstable gas particles in regions of converging flows. The star 
formation efficiency parameter e was set to 0.15, but with the 
adopted scheme the exact value of has only a minor effect on 
the star formation rate (Katz 1992). SN enrich the surround- 
ing gas according to the respective metal yields. Star particles 
inherit the metal abundance of the parent gas particles. How- 
ever we did not allow for diffusion of metals between gas par- 
ticles. 

After a gas particle is less than 10% of its initial mass due 
to multiple star formation events, it is removed and its mass 
is re-allocated among its gas neighbors. Up to six star parti- 
cles particles are then spawn for each gas particle in the disk. 
The adopted "minimal" feedback recipe dumps energy from 
SN into thermal energy of the nearest 32 gas neighbors. As 
energy is quickly radiated away in dense gaseous regions this 
feedback recipe has a weak effect on star formation (TC). We 
explicitly choose a fairly massive galaxy and of a weak feed- 
back to explore a regime where feedback is likely to play a 
minor role. 

The evolution and strength of the uniform UV background 



from QSOs followed Haardt and Madau (1996) and Haardt 
(2002, private communication). To complement our study, 
we performed dark matter only simulations of this same halo 
at similar and eight times higher resolution and 50% better 
spatial resolution. These simulations showed the main prop- 
erties of the galaxy's dark matter halo (inner profile, spin) are 
not changed substantially by the introduction of WDM, as the 
galaxy mass is much greater than the free-streaming mass at 
1O 95 M . 

3. RESULTS 

Within the ACDM galaxy virial radius (defined by Sp/p cr i t 
= 97.1 at z = 0) there are ~ 1 .2 x 10 5 dark matter, 5 x 10 4 gas 
and 6 x 10 5 star particles. The main parameters of the galaxies 
at z = are described in Table 1 . Star formation (Fig. 5) begins 
at z = 10, before the external UV background kicks in at z ~ 7. 
The star formation rate (SFR) peaks around z = 5 at 80MQ/yr 
while for the AWDM run the peak is delayed and is only 45 
M Q /yr at z = 4 due to the effect of removal of power at small 
scales. 

In both cosmologies the bulge of the galaxy forms through 
a series of rapid major merger events that end by z = 2.7, turn- 
ing mostly low angular momentum gas into stars. At z = 2.5 
the hot gas (T> 10 5 K) left within 50 kpc has a higher specific 
angular momentum than the dark matter (van den Bosch et al. 
2002). Soon after the last major mergers, the disk starts form- 
ing from the infall of this high angular momentum gas. At this 
stage the disk of the WDM galaxy is thinner than that formed 
in the ACDM one. In both simulations the final 3-4 Gyrs of 
evolution are marked by a bar instability, seen also in our own 
Milky Way (Cole & Weinberg 2002), and repeated encounters 
of the galaxy with a few relatively massive infalling satellites. 

3.1. The Galaxy Bulge and Disk components 

The resulting galaxies show that we successfully formed 
extended and fairly old stellar disks in A-dominated hierarchi- 
cal models. In Fig.l and Fig. 2 the stellar spheroid (including 
bulge and halo stars) is defined as the sum of all stars born 
before or at the time of the last major merger; this selects low 
angular momentum stars for which v rot j a < 0.5 (where v rot is 
the tangential velocity and a is the radial velocity dispersion), 
with the stars in the central few kpc having a v rot j a as low as 
0.1-0.2, similar to what inferred for the stars belonging to the 
spheroid of the Milky Way (Minniti 1996). We were not able 
to use a kinematic criterion in order to separate the central 
"bulge" from the diffuse extended spheroid — this resembles 
what found for our own Milky Way, where only a combina- 
tion of metallicities and kinematics allows a distinction (Min- 
niti 1996). We fit a de Vaucouleurs profile, characterized by a 
scale length Rb, to the old stellar distribution, and then we de- 
fine the mass of the bulge as twice the mass contained within 
Rb (Rb by definition corresponds to the half-mass radius of 
the system if we rely on the fit). This yields the total bulge 
masses indicated in Table 1 (precisely the ratio between the 
mass of the bulge and the mass of the disk, B/D is indicated 
there), and the bulge scale length is about 800 pc and 500 pc 
in, respectively, the ACDM and the AWDM run. 

At z = a flat, rotationally supported stellar component with 
Vmt I 'cr > 1 extends out to almost 20 kpc (Fig.l and Fig. 2). 
This disk component is thinner and dynamically colder (lower 
v w t /<?) in the AWDM model. The decomposition of the stel- 
lar distribution into a dynamically separate disk and a bulge is 
complicated by the presence of a third component, namely a 
prominent bar-like structure. The latter is significantly more 
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elongated and easier to distinguish in the AWDM galaxy. In 
the ACDM galaxy a "fat" bar or oval distortion typical of 
early-type galaxies (Athanassoula 2003) is first seen at z = 0.6 
and then is progressively thickened by tidal heating from in- 
falling satellites as well as by numerical heating, becoming 
progressively harder to distinguish. However, its signature is 
clearly evident in the kinematics; if we select all stars younger 
than 10 Gyr, thereby born after the last major merger, the 
kinematics are disk-like from 20 kpc going toward the cen- 
ter (v ' rot I a > 1.5, but become hotter (v ro ,/cr ~ 1) at 3-4 kpc 
from the center, the transition being sharper in the case of the 
AWDM galaxy (Fig.2). 

As mentioned above, the stellar density profiles of both 
galaxies at z = is thus the sum of the old spheroid formed by 
the last major merger and the much younger stellar bar in the 
inner few kpcs, whereas the cold disk dominates in the outer 
part (Fig.l and Fig.2). A double exponential provides a very 
good fit to the data. However, as from Fig. 1 it is worthwhile 
noting that a straightforward B/D decomposition based on the 
inner and outer component resulting from such fit would give 
erroneously large bulge masses and sizes, counting the thin, 
rotating and relatively young stars in the bar as part of the 
bulge. The overestimate of the bulge mass would be particu- 
larly severe in the AWDM galaxy. 

Given the initial halo spin (A = 0.03-0.035) the disks have 
scale lengths comparable to those predicted by theoretical 
models in which baryons conserve their specific angular mo- 
mentum during infall (see Mo, Mao and White 1998, which 
keeps into account the adiabatic contraction of the dark mat- 
ter halo due to the baryonic build up). The average age at the 
solar radius is 9-10 Gyrs for the ACDM run (8 for the AWDM 
run), in agreement with current observational estimates (e.g. 
van den Bergh 1996) In Fig.l we also plot the surface density 
profile of the old spheroid component - clearly this accounts 
only in part for the inner steepening. The bar contributes for 
the remaining part; the connection between the presence of a 
bar and a steep inner stellar profile has been widely shown by 
a number of numerical simulations on bar formation and evo- 
lution (Sparke & Sellwood 1987; Pfenniger & Friedli 1991, 
Mayer & Wadsley 2003). Interestingly, even in the case of the 
Milky Way the debate is still on-going on whether the central, 
kinematically hot component usually identified as the bulge 
is in part or even entirely a stellar bar (Zhao, Rich & Biello 
1996). 

The final spiral galaxies are more massive than M3 1 and the 
Milky Way. In both cosmologies the circular velocity (defined 
as v e = (M(r<R)/R)'/ 2 ) at the virial radius is ~ 185 km sec" 1 , 
while recent estimates (Klypin, Zhao & Somerville 2002) 
yield a v f at the virial radius around 150-160 km sec" 1 for our 
MW. At at 2.2 disk scale lengths (a distance typically used for 
observational samples) V2.2 = 320 km sec" 1 (280 km sec" 1 for 
AWDM, (see Fig. 3 and Fig. 4). As in other works, some of 
them including strong SN feedback (Abadi et al. 2002) the 
rotation velocity peaks even higher and in regions very close 
to the galaxy center that are only poorly resolved. These un- 
realistic high central velocities, linked to a very concentrated 
spheroidal component, are present in runs with and without 
SN feedback. They are a problem of most current galaxy for- 
mation runs. 

By the present time, star formation has turned about 64% 
(55% for AWDM) of gas inside the virial radius into stars, 
with ACDM having a slightly larger bulge component. At the 
final time only a small amount of gas within the virial radius 
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FIG. 3. — Rotation curve (denned as V c = ^fM({r < R)/R) for all compo- 
nents (continuous red line), star (blue) and dark matter (black) for the ACDM 
galaxy at z=0. The stellar kinematics are also shown; the rotational velocity 
is indicated by the green crosses, while the (azimuthally averaged) radial ve- 
locity dispersion in cylindrical coordinates (the z-axis being chosen along the 
angular momentum vector of the disk) is magenta 
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FIG. 4.— Same as for Fig.3, but for AWDM 



remains in the cold phase (1.5% compared to 2% at z=0.25 
and 9% at z=l in the ACDM). Most of it is in the central disk. 
The AWDM has three times as much cold gas in the disk. 
The exponential fits to the stellar disk mass distribution yield 
scale lengths of ~ 3 kpc for the ACDM and ~ 3.6 kpc for the 
AWDM. Bulge scale lengths are considerably smaller (Fig.l). 
The latter fits to the stellar distribution are obtained exclud- 
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ing the inner, bar-dominated region from the disk. Strictly 
speaking the bar is part of the disk since it was formed from 
disk material that underwent non-axisymmetric instabilities. 
Therefore we also tried fitting a single exponential over the 
whole stellar component with age < 10 Gyr; given the re- 
markable change of shape of the stellar density profile at a 
few kpcs from the center a single exponential provides a poor 
description. Nevertheless, it is worth noticing that, by doing 
this, the scale lengths drop to 2 kpc and 1.5 kpc for, respec- 
tively, the ACDM and the AWDM galaxy (as expected the 
change is more marked for the AWDM galaxy, which exhibits 
a stronger and longer bar). In the existing literature (see Car- 
ollo et al. 2001, Debattista,Carollo, Mayer & Moore 2003) 
the flatter component of the stellar profile is always included 
in the definition of the disk, while bars are often identified 
with a central bulge. These results support our double fits and 
the resulting scale lengths. 

By measuring the ratio of the specific angular momentum 
inside the virial radius of the halo in a dark matter only version 
of the simulations to that of the stellar disk at z = 0, we found 
that ~ 60% of the specific angular momentum has been lost in 
the ACDM galaxy and only 10% in the AWDM galaxy. The 
expected disk scale-length, R^, based on the models by Mo, 
Mao & White (1998) (see also Fall (1983) and Navarro & 
Steinmetz (2000)) is 3.4 kpc (ACDM) assuming a 50% loss 
of the specific angular momentum, thus very close to what 
we find, once the halo's spin and the adiabatic contraction of 
the dark halo due to the formation of both the bulge and the 
disk are taken into account. By neglecting the latter a much 
larger disk scale length, around 5 kpc, would be predicted 
by the analytical approach. A comparison between analytical 
predictions by Mo, Mao & White (1998) and observed disk 
galaxies (Courteau 1997) shows that the scatter in disk scale 
length measured by Courteau (1997) is more than a factor of 
3 across the entire circular velocity range. A galaxy with a 
scale length of about 3 kpc and V2.2 320 km/s stays lower than 
the mean on the relation predicted by the analytical model of 
Mo, Mao & White (1998), but still quite consistent with such 
data sample which consists mostly of Sb and Sc. However, 
to match the average properties of disk galaxies of type Sb 
and later a smaller angular momentum loss, or a higher initial 
halo angular momentum will be required. Indeed the ACDM 
galaxy would most likely be classified as Sa or Sab due to its 
prominent bulge. 

Using age and metallicity information from every star par- 
ticle created and Single Stellar Population models (Tantalo, 
Chiosi & Bressan 1998) we find that at z = the ACDM 
galaxy has M/ = -24.3 (including disk and bulge) and lies 
close to the Tully-Fisher relation (Giovanelli et al. 1997), 
Its disk B-R color is 1.05, while the average B-R in Courteau 
(1997) sample is 0.8. 

Internal instabilities of the stellar disk and perturbations 
by satellites cause structural and kinematical changes in the 
disk, heating it substantially after most of its mass has been 
assembled by z> 0.7-0.5. As mentioned above, the inner, 
colder regions where Q (Toomre's instability parameter) is 
initially < 2, become unstable to a bar-like mode at z ~ 0.6; 
part of the gas loses angular momentum to the stellar bar, 
flows to the center and becomes the site of a final episode 
of star formation. The disk of the AWDM galaxy stays thin- 
ner because tidal heating by satellites is largely suppressed. 
We tested (section 4) that the ACDM galaxy becomes bar- 
unstable mostly because of its fairly massive disk (Efstathiou, 
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FIG. 5.— SFR of the ACDM (continuous) and AWDM (dashed) galaxies. 
The line with small circles shows the SFR of the ACDM disk. 
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Lake & Negroponte 1982, Athanassoula & Sellwood 1986). 

3.2. The Galaxy Halo 

As expected in ACDM models, a rich amount of substruc- 
ture develops and forms a substantial stellar component due 
to a combination of weak feedback and early star formation. 
In excess of what is observed for Andromeda or the Milky 
Way (Mateo 1998) the ACDM model has at least thirty satel- 
lites with stellar masses larger than 1.35 xl0 9 M Q . Consistent 
with its lack of small scale power AWDM galaxy only has 
seven. 20% of the stars formed within the virial radius of the 
ACDM galaxy are in the satellites, and another 15% are in the 
halo within r>20 kpc. In comparison, due to its lack of power 
at small scales only 15% of stars are in the AWDM satellite 
system (3 of them are within 35 kpc of its center) and only 
0.2% of the stars are in the outer part of the halo. The most 
massive (and best resolved) satellites show a large variety of 
star formation histories (Fig. 6), in both cosmologies. Their 
star formation ceases a couple of Gyrs after entering the halo 
of the primary galaxy. At z = the stellar component of the 
dwarfs that entered the galaxy halo at z > 0.5 are supported 
by velocity dispersions, while those outside or that have just 
fallen in have recognizable disks, as expected if dynamical 
encounters are the major drivers of satellites morphologies as 
in the "tidal stirring" scenario (Mayer et al. 2001). In both 
models, the diffuse halos contain the oldest stars, but they 
also grow continuously from material stripped after the last 
major merger that gives rise to the bulge and from the debris 
of disrupted satellite until late epochs. 
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FIG. 6. — Star Formation histories for three satellites of the ACDM galaxy 
at z=(). Arrows mark the time they entered within R vi >. Satellites' SFRs peak 
in correspondence of merger events and pericentric passages 



The amount of baryons (including stars, cold and hot gas) 
inside the virial radius is within a few percent of the original 
cosmic abundance. The X-ray flux at z=0 is ~ 0.5 10 42 erg/sec 
in both cosmologies, this value doubles if a metal abundance 
of 0.3 solar is assumed for the gas in the halo. This is more 
than expected from observations, for example the large disk 
galaxy NGC2841 has a bolometric emission of only 1.3 10 41 
(Benson et al. 2000, Murali 2000). As expected with the 
weak feedback adopted and the long cooling times in the outer 
part of this relatively massive halo a hot gas phase is always 
present within the virial radius, and X-ray luminosity evolu- 
tion is weak to z=l. The overabundance of halo baryons (in 
all its various subcomponents, hot gas, satellites and diffuse 
stellar halo) is clearly a problem of both cosmologies. While 
the WDM galaxy has a substantially smaller stellar halo com- 
ponent, possibly in agreement with observations, its X-ray 
emission is high and very similar to theACDM one. This is 
expected as reduced power at subgalactic scales and does not 
affect the total mass and the accretion rate of halo baryons. 
According to the existing literature and the work presented 
here, energy feedback, reduced amount of small scale power 
and increased resolution tend to increase the size of galac- 
tic disks, but the amount of hot gas in halos of large galax- 
ies is likely only sensitive to feedback. X-ray observations 
sensitive to sub keV temperatures of nearby massive galax- 
ies sensitive would set powerful constraints on the amount of 
energy injected by different processes into the baryonic com- 
ponent during the formation of the galaxy, as they currently 
do on the larger scale of small groups (Borgani et al. 2002), 
pointing to at least 0.5keV/baryon. As an alternative Toft et 
al. (2002) suggest that the excess X-ray emission could be re- 
duced allowing a detailed treatment of cooling by metal lines 
that would reduce the amount of hot gas left in the halo. Cur- 
rent global estimates of galaxy formation efficiency set it at 
about 10% (Bell et al. 2003) which is significantly lower than 



the high star formation efficiency in our runs, suggesting that 
ejection of baryons from the halos of disk galaxies is likely 
necessary to avoid excess X-ray emission from galaxy halos. 
The hot gas component of the Milky Way has just recently 
been detected (see Sembach 2003) and future observations 
will provide useful constraints on its properties. 

4. ROBUSTNESS OF NUMERICAL RESULTS: SOLVING THE 
ANGULAR MOMENTUM PROBLEM 

The main finding of this work is that the angular momentum 
and consequently the scale length of the stellar disk formed in 
our simulations is comparable to that of real galaxies of simi- 
lar luminosity. This result was not obtained adopting the usu- 
ally advocated strong feedback to decrease the lumpiness of 
the gas being accreted, but rather using good mass and force 
resolution. Is then the addition of feedback still a necessary 
ingredient to obtain realistic galaxies? How much is the for- 
mation and evolution of disk stellar systems affected by lim- 
ited resolution? 

It is now widely recognized that 10 5 " 6 particles with a force 
resolution close or better than 1% R„> are necessary to re- 
liably simulate the formation and the internal structure of 
a single dark matter halo (Ghigna et al 2000, Fukushige & 
Makino). It is then unlikely that the much more complex 
formation of disk galaxies could be reliably simulated with 
significantly less particles. A comparison with our higher res- 
olution dark matter only runs shows that the mass function 
of dark matter subhalos becomes severely incomplete below 
30-35 km sec" 1 . On the positive side this implies that at our 
resolution we are resolving the progenitors mass range that 
contributes most of the mass in the accretion history of the 
central galaxy and most of the mass in the subhalo popula- 
tion. With more than 10 5 particles within the virial radius, our 
simulations accurately follow the evolution of the hot gaseous 
phase (see tests in Frenk et al. 1999, Borgani et al. 2002, 
Wadsley, Stadel & Quinn 2003). 

The final angular momentum of a stellar galactic disk is set 
first by the details of hierarchical accretion onto the parent 
halo and the subsequent cooling of high angular momentum 
gas. Once a disk of baryonic material is formed, disk instabil- 
ities might redistribute the original angular momentum of the 
disk both internally and to the dark matter halo (Weinberg & 
Katz, Valenzuela & Klypin 2003). Both stages are potentially 
affected by numerical viscosity (present in the SPH approach 
to describe shocks) and numerical relaxation introduced by 
noise in the global potential due to small particle numbers. 
These effects reduce the angular momentum of the disk by 
artificial dissipation or by transfer to the halo component re- 
spectively. 

To measure the importance of the numerical effects men- 
tioned in the above paragraph, but independently of the hier- 
archical build up, we built an isolated galaxy system (follow- 
ing Springel & White 1999) with the same structural param- 
eters (mass, profile, fraction of mass in the disk and in the 
bulge, relevant scale lengths, fraction of cold gas in the disk) 
as our ACDM galaxy at z=0.6, i.e. soon after the formation 
of the disk. Results can likely be extended to the AWDM 
galaxy. All components were made of "live" particles This 
galaxy was evolved in isolation for 6 Gyrs. We run different 
version of it changing the number and accordingly the mass 
and softening of the stellar, dark matter and gas particles (see 
Fig. 7). Softening for the different component in the different 
realizations was scaled as e test = e cosmo (m^sf/m,^^) 1 / 3 , were 
(■cosmo and m cosm „ were the softening and masses used in the 
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FIG. 7. — Angular momentum loss in an isolated disk galaxy model with 
structural parameters as our ACDM run at z = 0.6. This model was run for 
6 Gyrs ( equivalent to the present time). The y axis shows the fractional 
angular momentum loss for all the baryonic material in the disk as a function 
of radius. Continuous line:N DM = 100000, Nstar=200000, Ngas=50()0 (same 
as in the cosmo run). Our choice of spawning a large number of star particles 
(typically 6 per gas particle) in the cosmological run reduces internal heating 
of the stellar disk, dotted short dashed: N/jm and Nstar reduced by a factor 
of five. Long dashed: Nstars reduced by a factor of 25. short dashed:NoM= 
4000, Nstar=8()00, Ngas=l()0(). At the lowest resolution the disk undergoes 
the catastrophic angular momentum loss reported in early simulations 



main cosmological runs for each component. 

In our lowest resolution realization (25 times less particles 
than in our cosmological runs), the number of halo (DM) par- 
ticles and particles representing stars in the disk goes below a 
few thousands and the disk (defined as all the particles origi- 
nally assigned to the disk component) loses almost 70% of its 
initial angular momentum over 8 Gyrs. A disk structure is not 
visible at the final (present) time (Fig. 7 short dashed). Angu- 
lar momentum loss is greatly reduced as the number of DM 
particles first and then number of particles in the stellar disk 
is increased (long dashed and dot dashed). At the resolution 
of our cosmological runs (continuous line) the loss is down to 
only 10% across most of the disk, reaching 30% in the inner 
2 kpc. Even in this isolated model a bar forms in the disk. 
Finally, we verified that when the collisionless components 
are well resolved then gas angular momentum loss due to ar- 
tificial viscosity becomes less than 10% over several Gyrs if 
the disk gas component is represented by at least a few thou- 
sand particles (in agreement with Navarro & Steinmetz 2000). 
We expect numerical effects will continue to decrease as the 
numerical resolution (both spatial and in time) is increased. 

To verify that our findings are applicable in a full cosmo- 
logical context where the hierarchical build is included, we 
also run a low resolution realization (l/8th of the particles 
and twice as large e of the cosmological ACDM galaxy de- 
scribed in our work. In this situation additional resolution 
problems could be present, as the first individual structures 
that form are necessarily made of few particles. Artificial vis- 



cosity could show up not only in cold gas lumps, but also in 
poorly resolved shocks and numerical effects could cascade 
down as the galaxy main body is assembled. This also im- 
plies that, at the same mass and force resolution, the ACDM 
galaxy with its larger substructure mass range, is intrinsically 
harder to simulate than the AWDM one. As expected, due 
to the low number of particles and the very poor spatial res- 
olution by z=0 in this low resolution run, cold gas and stars 
within 20 kpc have lost 90% of their initial angular momen- 
tum compared to the dark matter and no strong stellar disk 
component is present, although the central part of the stellar 
distribution is rotationally supported. This low resolution run 
is in qualitative agreement with Eke, Navarro & Steinmetz 
2000, where no disk like morphology was observed in their 
simulated galaxies. Our results imply that the disk galaxy de- 
scribed in TC01 (their run was stopped at z=0.5) would likely 
show a much weaker disk component if evolved to the present 
time. 

To show the effects of resolution on a galaxy with differ- 
ent mass and merging history from that in our main runs we 
re-simulated to z= 2 (with similar techniques as described for 
the main runs) a small dwarf galaxy sized halo (9x 1O 1O M ) 
at the outskirts of our MW halo at 8 times better mass reso- 
lution. As the number of DM particles within R vlr (gas parti- 
cles were increased accordingly) increased from ~ 3000 to 
almost 2.5 10 4 and the softening was set at 0.5 kpc, the 
stellar component of this dwarf galaxy became more rota- 
tionally supported as v/er changed from 0.46 to 0.77 and a 
gaseous disk became apparent. The shape and amplitude of 
the SFR(z) did not change substantially, showing that we are 
qualitatively capturing the events that drive star formation in 
the small galactic satellites of our main runs. 

To summarize, at the best resolution for our main runs nu- 
merical effects are likely responsible for the angular momen- 
tum loss in the AWDM disk (about 10%). As described in 
section 3 the ACDM disk loses significantly more angular 
momentum than its equivalent in AWDM. The only differ- 
ence in the two runs are the initial power spectra, with ACDM 
having significantly more substructure below 50-100Km/sec. 
This result suggests that at the current resolution, and as sug- 
gested by early works the intrinsic dumpiness of the hier- 
archical formation process leads to significant angular mo- 
mentum loss, likely through interactions between the infalling 
lumps (Lake & Carlberg 1986) before the disk has formed. 
However, at our resolution, this is angular momentum loss 
is far from catastrophic as the fraction of mass contained in 
surviving sublumps within a larger parent halo is only of the 
order of 15% (Ghigna et al 1998). A large fraction of mass is 
tidally stripped from the subhalos as soon as they first go past 
pericenter (Mayer 2001). This material typically loses only a 
small fraction of its initial angular momentum (Colpi, Mayer 
& Governato 1999.) If this high fraction of stripped material 
is also typical of the gas accreted by the disk then most of 
it actually forms from diffuse material cooling from the hot 
halo, rather than by individual lumps. 

Our results suggest that the large angular momentum loss 
or the lack of well defined stellar disks reported in early 
works might have resulted from poor force or mass resolu- 
tion. The tests done with the isolated galaxy model suggest 
that two-body heating by massive halo particles would be the 
main cause of angular momentum loss at low resolution if 
small softenings were used. As done in some previous works, 
two-body heating can be suppressed, by using a larger soft- 
ening without increasing the number of particles. However, 
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even softenings as small as those values used here are barely 
enough to resolve the disk and avoid spurious results as for 
its morphological evolution; suppression of disk self-gravity 
by a softening larger than the disk scale length would inhibit 
processes like bar formation that are essential ingredients in 
disk formation and evolution (Romeo 1991, 1994; Weinberg 
& Katz 2002; Valenzuela & Klypin 2003; Mayer & Wadsley 
2004). 

We also evolved another isolated galaxy model with the 
same structure and numerical resolution as the ACDM galaxy, 
but with a stellar disk only 50% as massive. This disk does 
not show the bar instability present in the other runs, and its 
angular momentum loss is slightly lower. This experiment 
also implies that disk-satellite interactions have only a lim- 
ited effect on the development of the bar instability and in the 
angular momentum loss of the disk. 

Finally, it is likely that our adopted feedback recipe is too 
simplistic, and although it models the build up of the disk pop- 
ulation well, it does not stop star formation at high redshift 
in small halos creating an excess of massive satellites in the 
ACDM galaxy The paucity of cold gas in our final state sug- 
gests that our star formation algorithm is likely too efficient, 
at least compared to the universal value of 10% or less (Bell 
et al. 2003) although is difficult to quantify such a statement 
with only one run. 

5. SUMMARY AND DISCUSSION 

We simulated the formation of a disk galaxy in a cosmo- 
logical context to the present time with a high force and mass 
resolution. Without resorting to strong feedback recipes, we 
successfully formed an object with bulge-to-disk mass ratios, 
scale lengths, current star formation, disk age and dynami- 
cal properties representative of those observed in present day 
massive spirals. Our results show that the large angular mo- 
mentum loss reported in some early works was in at least in 
part due to insufficient mass or force resolution. High resolu- 
tion plays a significant role in simulating correctly the forma- 
tion of disk galaxies in a cosmological context, independently 
of the SN feedback recipe adopted. 

We estimate that at a resolution of ~ 10 5 particles within 
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